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ABSTRACT 

A sample of 47 faint Gigahertz Peaked Spectrum (GPS) radio sources selected from the 
Westerbork Northern Sky Survey (WENSS, Rengelink et al. 1997), has been imaged 
in the optical and near infrared, resulting in an identification fraction of 87%. The 
R—I and R—K colours of the faint optical counterparts are as expected for passively 
evolving elliptical galaxies, assuming that they follow the R band Hubble diagram 
as determined for radio-bright GPS galaxies. We have found evidence that the radio 
spectral properties of the GPS quasars are different from those of GPS galaxies: The 
observed distribution of radio spectral peak frequencies for GPS sources optically 
identified with bright stellar objects (presumably quasars) is shifted compared with 
GPS sources identified with faint or extended optical objects (presumably galaxies), 
in the sense that a GPS quasar is likely to have a higher peak frequency than a 
GPS galaxy. This means that the true peak frequency distribution is different for 
the GPS galaxies and quasars, because the sample selection effects are independent 
of optical identification. The correlation between peak frequency and redshift as has 
been suggested for bright sources has not been found in this sample; no correlation 
exists between R magnitude (and therefore redshift) and peak frequency for the GPS 
galaxies. We therefore believe that the claimed correlation is actually caused by the 
dependence of the peak frequency on optical host, because the GPS galaxies are in 
general at lower redshifts than the quasars. The difference in the peak frequency 
distributions of the GPS galaxies and quasars is further evidence against the hypothesis 
that they form a single class of object. 



1 INTRODUCTION 



Extragalactic radio sources have been studied for many 
years, but it is still unclear how they are formed and evolve. 
A crucial element in the study of their evolution is to iden- 
tify the young counterparts of the "old" FRI/FRII extended 
objects. Good candidates for young radio sources are those 
with peaked spectra (Gigahertz Peaked Spectrum (GPS) 
sources and Compact Steep Spectrum (CSS) sources, Fanti 
et al. 1995, Readhead et al. 1996, O'Dea and Baum 1997, 
O'Dea 1998), because they are small in angular size as ex- 
pected for young sources. GPS sources are characterized by 
a convex radio spectrum peaking at a frequency of about 1 
GHz and are typically 100 pc in size. CSS sources have a 
peak in their spectrum at lower frequencies and have pro- 
jected linear sizes of < 10 — 15 Kpc. A commonly discussed 
alternative to them being young, is that they are small due 
to confinement by a particularly dense and clumpy inter- 



stellar medium that impedes the outward propagation of 
the jets (O'Dea et al 1991). 

Important information about the nature of Gigahertz 
Peaked Spectrum GPS sources comes from the properties of 
their optical counterparts. Systematic optical identification 
programs on GPS sources have been performed by O'Dea 
et al. (1990), Stanghellini et al. (1993) and de Vries et al. 
(1995), and these show that the optical counterparts are a 
mixture of galaxies and quasars. O'Dea et al. (1990) and 
Stanghellini et al. (1993) found that disturbed optical mor- 
phologies are a common characteristic of GPS galaxies, with 
many of them having close companions or appearing to lie 
in a group or cluster. This may indicate that galaxy-galaxy 
interactions and/or mergers play an important role in the 
GPS phenomenon. GPS galaxies appear to be a homoge- 
neous class of giant ellipticals with old stellar populations 
(Snellen et al. 1996a, 1996b, O'Dea et al. 1996). Their opti- 
cal to near-infrared fluxes seem not to be contaminated by 
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light produced by the active nucleus, and they are therefore 
potentially useful probes of galaxy evolution. 

In the light of the orientation unification scheme (eg. 
Barthel 1989), it would be expected that the redshift distri- 
butions of GPS sources identified with galaxies and quasars 
are more or less the same. However, the GPS quasars tend 
to be at much higher redshift than the galaxies (2 < z < 4, 
O'Dea et al. 1991). In addition, the intrinsic peak frequency 
distributions and radio morphologies seem to be different for 
GPS galaxies and quasars (Stanghellini et al. 1996). These 
results are probably influenced by several selection effects, 
but if true, they may indicate that GPS galaxies and quasars 
are not unified by orientation, but form separate classes of 
object which happen to have similar radio spectra. 

To further address these issues, we have started to in- 
vestigate GPS sources at fainter flux density levels (Snellen 
et al. 1995, Snellen et al. 1998) than those in the samples of 
bright GPS sources studied until now (O'Dea et al. 1991, de 
Vries et al. 1997, Stanghellini et al. 1998) . By comparing the 
properties of these new faint samples with those of bright 
samples in the literature, we shall investigate properties of 
GPS sources as a function of radio luminosity, redshift and 
rest frame peak frequency. 

This paper describes the optical and near-infrared imag- 
ing of a sample of faint GPS sources selected from the 
Westerbork Northern Sky Survey (WENSS Rengelink et al. 
1997). 



2 THE SAMPLE OF FAINT GPS SOURCES 

The selection of the sample has been described in detail in 
Snellen et al. 1998 (paper I). Candidate GPS sources selected 
from the Westerbork Northern Sky survey, are those with an 
inverted spectrum between 325 MHz and higher frequencies. 
The sources are located in two regions of the survey; one at 
15 h < a < 20 h and 58° < 5 < 75°, which is called the 
mini-survey region (Rengelink et al. 1997), and the other at 
4 h 00 m < a < 8 h 30 m and 58° < 5 < 75°. Additional ob- 
servations at 1.4, 5, 8.4 and 15 GHz were carried out with 
the Westerbork Synthesis Radio Telescope (WSRT) and the 
Very Large Array (VLA), yielding a sample of 47 genuine 
GPS sources with peak frequencies ranging from 500 MHz 
to more than 15 GHz, and peak flux densities ranging from 
~ 30 to ~ 900 mjy. To determine whether the GPS sources 
have an optical counterpart on the Palomar Sky Survey 
plates (POSS), the Automatic Plate Measuring (APM) cat- 
alogue (Irwin et al. 1994) and the HST Guide Star Catalogue 
(Russel et al. 1990) were used for the sources in the mini- 
survey region and in the region 4 h 00 m < R.A. < 8 h 30 m , 
respectively. This information was used for finding charts, 
to determine the required exposure times for the observa- 
tions, and for astrometric purposes. 



3 OBSERVATIONS 

The sample was imaged in the optical in the course of 
three observing sessions (see table [j]), with the 0.9m Jacobus 
Kapteyn Telescope (JKT), the 2.5m Isaac Newton Telescope 
(INT) and the 2.5m Nordic Optical Telescope (NOT), all 
located at Observatorio del Roque de los Muchachos on La 



Palma, Spain. The sample was also observed in the near- 
infrared with the 4.2m William Herschel Telescope (WHT), 
and with the NOT. The INT, NOT and WHT observations 
were done as part of the international CCI observing time 
awarded to the WENSS team in 1995. 

During the different observing sessions, several standard 
stars from Landolt (1983) were observed for photometric 
calibration. A summary of the relevant parameters for the 
observations is given in table 

JKT CCD Imaging 

Six nights of observing time were obtained on the JKT, 
from 3 to 9 June 1994. All 25 GPS sources in the sample 
from the mini-survey region were observed in R and I bands 
using Harris filters (Matura et al. 1993), which have, in com- 
bination with the GEC7 CCD, effective wavelengths of 6550 
and 8260 A and FWHM bandwidths of 1500 and 3000 A for 
the R and I filters respectively. Typical observing times were 
3 x 300 seconds for the POSS detections and 3 x 600 sec- 
onds for the unidentified sources on the POSS. The 400 x 590 
GEV7-chip used as the detector has a projected pixel size 
corresponding to 0.3 ", and a total field of view of 2' x 3'. In 
addition to the Landolt standards, spectroscopic standards 
were measured to determine the absolute flux density scale 
for the particular set of filters. Limiting surface brightnesses 
(lcr) were for the longest exposures typically 26.1 and 25.3 
mag/arcsec -2 in respectively R and I band. The weather 
was photometric during the run, except for the second half 
of the night on 7 June. The sources observed on 7 June were 
re-observed briefly on 9 June for photometric calibration 
purposes. 

INT CCD Imaging 

Sources which had not been detected with the 
JKT were reobserved using the INT on 6 April 1995. 
These objects, B1525+6801, B1551+6822, B1600+7131, 
B1620+6406, B1639+6711, B1655+6446, and B1808+6813, 
were observed in R band using a 1024x1024 TEK-chip and 
Harris filter (A e// = 6550A, FWHM=1500 A). Each source 
was observed for 15 minutes under photometric conditions. 
Limiting brightnesses (lcr) were typically 27.3 mag/arcsec -2 
for the longest exposures. The CCD detector had a pixel 
scale of 0.58" and a field of view of 10' x 10'. 

NOT CCD Imaging 

Sources in the region 4 h 00 m < R.A. < 8 h 30 m and 58° < 
decl. < 75° were observed using the NOT in the period 1 
to 5 October 1995. A 1024x1024 TEK Chip in combination 
with a Harris R band filter (A e// = 6550A, FWHM=1500 
A) was used. The CCD detector had a pixel scale of 0.17" 
and a field of view of 3' x 3'. One source, B1954+6146, was 
also observed using an I band filter which, in combination 
with the TEK CCD, had an effective wavelength of 8300 A. 
Limiting brightnesses (la) were typically 27.0 mag/arcsec -2 
for the longest exposures. The conditions were photometric 
during the whole run. 

NOT ARNICA Observations 

About half of the GPS sources in the mini-survey region 
were observed in K band (A = 2.20 ± 0.20/im), and some in 
J band (A = 1.28±0.15/im)using the ARcetri Near Infrared 
CAmera (ARNICA) on the NOT, from 15 to 17 August 
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Telescope 


Observing 


Detector 


Filter 


Detector 


Pixel 


Field of 




Date 






Size 


Size 


View 


JKT 


3-9 Jun. 1994 


CCD 


R, I 


400 x 590 


0.30" 


2' X 3' 


INT 


6 Apr. 1995 


CCD 


R 


1024 X 1024 


0.58" 


10' x 10' 


NOT 


1-5 Oct. 1995 


CCD 


R 


1024 X 1024 


0.17" 


3' x 3' 


WHT 


28&31 Jan. 1996 


WHIRCAM 


K 


256 x 256 


0.24" 


1' x 1' 


NOT 


15-17 Aug. 1995 


ARNICA 


J, K 


256 x 256 


0.29" 


1.25' x 1.25' 


NOT 


5-6 Sep. 1995 


ARNICA 


J, K 


256 x 256 


0.29" 


1.25' x 1.25' 



Table 1. Details of the optical and near- infrared observations. 



and on 5 and 6 September 1995. These observations are 
discussed in detail by Villani and di Serego Alighieri (1998). 
ARNICA utilizes a NICMOS 3 (256 x 256 pixels) array 
as a detector, with a projected pixel size corresponding to 
0.29" and a field of view of 75" x 75". In order to ensure 
background limited images and to avoid saturation of the 
brightest sources, the total exposure time on each source was 
divided into a number of shorter exposures, between which 
the telescope was moved in a 9 position raster. The total 
exposure times were between 6 and 30 minutes, depending 
on the expected magnitude of the object. 

WHIRCAM Observations 

Seventeen sources in the region 4 h 00 m < a < 8 h 30 m 
were observed on 28 and 31 January 1996, using the WHT 
Infrared Camera (WHIRCAM). WHIRCAM utilizes a 256 x 
256 InSb array detector, and was used as a direct camera 
with an image scale of 0.24 arcsec/pixel. Each source was 
observed in K short band (A = 2.16 ± 0.16/im) with an 8 
position jitter pattern. At each position the field was inte- 
grated for 2, 4 or 10 x 10 seconds using the ND.STARE 
mode, meaning that the array is reset and read immedi- 
ately and read again after the exposure time. This resulted 
in total exposure times of 160, 320 or 800 seconds, depend- 
ing on the expected magnitude of the object. Two sources, 
B0830+5813 and B0755+6354, were observed on 28 January 
1996 through thin cirrus. The conditions on 31 January were 
photometric. 

3.1 Optical Data Reduction and Calibration 

The reduction of the optical observations was carried out 
with the software packages IRAF from the NOAO and AIPS 
from the NRAO. The fits images were loaded into IRAF and 
bias subtracted and flat fielded in the standard manner. Po- 
sitional calibration was then carried out within AIPS using 
the task XTRAN. Positions of stars in the CCD fields were 
taken from the APM catalogue (Irwin et al. 1994) for the 
sources in the mini-survey region and the HST Guide Star 
Catalogue for the sources with 4 h 00 m < a < 8 h 30 m . These 
positions have an intrinsic rms accuracy of typically 0.7". 
The epoch and equinox of the reference coordinate system 
was then converted from B1950 to J2000 using the AIPS 
task EPOSWITCH. 

The absolute flux density scale was determined using 
the spectroscopic standard stars observed during the JKT 
run. This was done by combining the tabulated stellar spec- 
tra (Stone 1977) with the filter relative transmission curves 
and the CCD response. An object of zero magnitude, ob- 
served in the Harris R band filter, was estimated to have 



a flux density of 2.2 ± 0.2 x 10~ 9 erg cm' 2 sec' 1 A' 1 at 
6550 A. An object of zero magnitude, observed in the Har- 
ris I band filter, has been estimated to have a flux density 
of 1.0 ± 0.1 x 10~ 9 erg cm' 2 sec' 1 A' 1 at 8900 A. The 
filter response curves indicate that these absolute calibra- 
tion measurements can also be applied to the Harris R band 
observations on the NOT and the INT, due to the similar 
CCD characteristics. The I band filter used on the NOT had 
no absolute photometric calibration information available. 
However, since the filter bandpass is similar to the Harris I 
band filter used on the JKT, the same absolute calibration 
was used. 

3.2 Near-infrared Data Reduction and 
Calibration 

The data of the ARNICA and WHIRCAM observations were 
reduced using IRAF. The images were dark-subtracted and 
then flat-fielded by median filtering the image stack. There- 
after the frames were normalized and mosaiced together to 
produce one final image per object. In the case of the AR- 
NICA observations, it was necessary to do the background 
subtraction from each frame before flat-fielding, in order to 
eliminate the contribution of the thermal emission from the 
background and the telescope. UKIRT faint standards were 
observed several times during the nights at a range of air- 
masses. Although the WHIRCAM observations were done 
using a Jfshort filter, the K band magnitude system has 
been used for which zero magnitude corresponds to a flux 
density of 4.1 ± 0.4 x 10~ n erg cm' 2 sec' 1 A' 1 . Astro- 
metric measurements have been done in the same way as in 
the optical, using positions of objects in the R band CCD 
images that were also in the near-infrared fields. 

3.3 Identification Procedure 

The procedure for identifying a GPS source in the optical is 
simplified compared to large size radio sources, because the 
angular size of the radio source is very small compared to the 
optical resolution. The likelihood method (eg. de Ruiter et 
al. 1977) was used to quantify our identification procedure. 
The reliability of a proposed identification has been assessed 
using the observed differences between the optical and ra- 
dio positions and their errors and the optical background 
counts. For each candidate counterpart the following dimen- 
sionless measure of the uncertainty in position difference, r, 
was calculated. 

r = J-s— + -s— (!) 

V a aopt a Sopt 
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where Aa and AS are the offsets between the optical and 
radio positions, and a aop t and asopt the uncertainties in the 
optical right ascension and declination positions. The error 
in the radio position is small compared to the error in optical 
position (0.2" versus 0.7") and was therefore neglected. 

An optical object is considered to be a candidate identi- 
fication if it is located within an angular distance 3.5" from 
the radio position. Only six sources did not have a candi- 
date identification in the optical or near-infrared. The radio- 
optical offsets are shown in figure |l| for the sources in the 
mini-survey region and in the region 4 h 00 m < a < 8 h 30 m . 
For the sources in the region 4 h 00 m < a < 8 h 30 m and 
58° < S < 75°, a median position offset of 0.6" in declina- 
tion was found between the optical and the radio. This was 
removed before determining the dimensionless position un- 
certainties. The cause of this offset is not clear. The median 
position offsets in a for the region 4 h 00 m < a < 8 h 30 m , and 
in both a and 5 for the sources in the mini-survey region 
were found to be smaller than 0.05" and therefore not taken 
into account. 

The background density of objects, pt g , on the CCD 
frames, depends mainly on the galactic latitude and mag- 
nitude limit of the observation, and was in each case esti- 
mated by counting the number of objects on the CCD frames 
for each source. The background densities range from 5 to 
~ 40 objects per square arcminute. The likelihood method 
assumes that a radio source and its optical identification al- 
ways coincide in position, and that the true identification is 
always closer to the radio source than the first contaminat- 
ing object. Given the normalized position difference r for a 
certain radio optical pair, the likelihood ratio is defined by 
(de Ruiter et al. 1977); 



identified sources for which the identification is correct, are 
given by: 



LR{r) 



dp(r\id) 



r e 



2A r e 



-Xr 2 



1 ,- 2 (2A-l) 

— e 2 
2A 



(2) 



dp(r\c) 

where A = TV(T a0 ptO'soptPbg = l-54pf, 9 . This is the ratio of the 
probability that a given object found between r and r + dr 
is the correct identification p(r\id), divided by the probabil- 
ity that it is a contaminating object p(r\c). We have used a 
likelihood ratio cut off of 2.0, which means that a source is 
considered to be identified if the probability that the given 
object is the correct identification is twice the probability 
that it is a contaminating background source. A likelihood 
ratio of 2.0 corresponds to an optical to radio offset of 2.9" 
for a background density of 10 objects per square arcminute. 
Note that all the candidate identifications fall within this 
limit, except for B1647+6225 which has two candidate coun- 
terparts, of which the closest one is selected as a possible 
identification. 

Using a likelihood ratio cutoff of 2.0, we find an iden- 
tification fraction 8 of 87%. This identification fraction can 
be related to the two a posteriori probabilities that the ob- 
ject found at an angular distance r from the radio source 
position is a genuine identification, p(id\r), or a confusing 
object, p(c\r), using Bayes' theorem (eg. Rottgering et al. 
1995). 

6 LR(r) 



p(id\r) — 



p(c\r) = 



(3) 



6 LR{r) +1-6* v 1 ' 9 LR{r) + 1 

The completeness, C, of the identifications, which is 
the number of accepted identifications over the number of 
correct id's, and the reliability, R, which is the fraction of 



C=l-( ^2 P{id\r))/N id 



R 



LR( ri )<L 

E 

LR(ri)>L 



P(c\r) )/N id 



(4) 



(•») 



Sources in our sample which remain unidentified have 
no optical and/or near-infrared counterpart within 3.5" of 
the radio position. This corresponds to a completeness of 
C > 99.9%. The reliability is determined to be 97.6%, which 
implies that it is probable that one source will have an er- 
roneous identification. 



3.4 Parameters of the Optical Counterparts 

The magnitudes enclosed within a circular aperture were 
measured for all the candidate identifications within var- 
ious radii from the peak brightness of the object. These 
radii are between 3.3" and 24", depending on the angu- 
lar size of the optical identification. Furthermore, for the 
sources with tur < 21.5, the FWHM of the object was mea- 
sured and compared to that of stars in the frame, in or- 
der determine if it is an unresolved or an extended source. 
A 3cr detection limit in a 3" aperture is taken as an up- 
per limit for sources without a candidate identification. The 
measured magnitudes were then corrected for the galactic 
foreground extinction using data from Burnstein and Heiles 
(1984) derived from HI column densities, assuming an ex- 
tinction curve proportional to A" 1 . This might be too large 
a correction in the near-infrared, however extinction correc- 
tion are small at those wavelengths anyway. These extinction 
data have been taken from the NASA/IPAC extragalactic 
database (NED). 



4 RESULTS AND DISCUSSION 

The results for the individual sources are presented in table 
^| Column 1 gives the source name; columns 2 the optical 
positions (J2000) for the candidate identifications; columns 
3 and 4 the radio-optical position offsets in right ascension 
and declination; column 5 the dimensionless position error; 
column 6 gives the morphological classification of the identi- 
fication (F=Faint R > 21.5, S=Star-like, E=Extended), col- 
umn 7 the galactic foreground extinction in B band, column 
8 the aperture diameter used for measuring the magnitudes 
and column 9, 10, 11 and 12 give the R, I, J, and K band 
magnitudes respectively. 

The optical and near-infrared images are presented in 
the appendix, with notes on individual sources. All the R 
band and WHIRCAM K band, and the I band images 
with detections are shown. The near-infrared maps taken 
with ARNICA are shown by Villani and di Serego Alighieri 
(1998). The R and / band images show a region of 1' x 1', 
and the K band images a region of 30" x 30" around the 
radio position. To reduce the effect of digitization of the 
contours, the data was convolved with a Gaussian with a 
FWHM of 3 pixels. The contours are at 3, 6, 12, 24, 48, etc. 
times the sky noise for the bright sources, and at 3, 4, 6, 12, 
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Figure 1. The radio-optical position offsets in arcseconds for the candidate identifications in the mini-survey region (left) and in the 
region A h 00 m < a < 8 h 30 m (right). The squares indicate candidate identifications for which there are other objects closer to the radio 
position. Note that there is a median offset in declination of 0.6" for the sources in the region 

4 h 00 m < a < g h 30 m , which has been 
removed before determining the dimensionless position errors. The dotted circles indicate the uncertainties of the optical position. 




14 16 18 20 22 24 26 

R-Band Magnitude 



Figure 2. The R band magnitude distribution of the sources in 
the sample. The shaded areas indicate identifications with star- 
like objects, which are assumed to be quasars. The arrows indicate 
lower limits. 

24, 48 etc. times the sky noise for the faintest sources. The 
level of the faintest contour in mag/arcsec 2 is given at the 
top-right corner of each image. 

4.1 Magnitude and Colour Distributions 

Figure ^| shows the R band magnitude distribution of the 
sources in the sample. If a radio source has more than one 
candidate identification, the one closest to the radio position 
was taken as the ID. The shaded areas indicate identifica- 
tions with star-like objects, which are assumed to be quasars. 
In our sample there are 6 bright galaxies (13±5%), 20 bright 
quasars (43±9%), and 21 faint objects (44±10%) for which 



morphological classification is uncertain. This was compared 
with the radio bright sample of Stanghellini et al. (1998) con- 
taining 33 sources, which is complete to a flux density of 1 
Jy at 5 GHz. For this radio bright sample, the percentages 
are 36±10% for the bright galaxies, 42±11% for the quasars, 
and 21 ± 8% for the faint objects with thr > 21.5. While 
the percentages of bright quasars in both samples are the 
same within the errors, the ratio of bright to faint galaxies 
is lower in our faint sample than in the radio bright sample 
at a 2a level. Assuming that the faint optical objects are 
galaxies, this may indicate that the redshift distribution of 
GPS galaxies in our sample is more biased towards higher 
redshift. 

Figure ^ shows the available R — I and R — K colours 
as function of R band magnitude. The diamond symbols 
represent the extended objects, the crosses the star-like ob- 
jects, and the squares the faint ma > 21.5 identifications. 
The dotted lines show the upper and lower limits on colour 
for a spectral energy distribution of F v on v ~ clo v t with 
0.0 > a opt > —2.0, which is the expected range for a opt of 
radio loud quasars (eg. Baker and Hunstead 1995). All the 
colours of the stellar objects are within the expected range 
for quasars, except for two R — I colours and one R — K 
colour. These are slightly redder than expected which may 
be due to the contribution of emission lines in the / or K 
band. 

The solid lines represent the R— I and R — K colours ex- 
pected for passively evolving ellipticals, assuming that the 
galaxy magnitudes lie on the R band Hubble diagram for 
GPS sources defined in Snellen et al. (1996a). The R - I 
colours of the extended and faint identifications are con- 
sistent with this, but four from the six faint objects with 
measured R — K colours are bluer (> 2a) than expected. 
These could still be passively evolving ellipticals with high 
formation redshifts as long as they are intrinsically fainter 
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Figure 3. R — I (top) and R — K (bottom) colours as func- 
tion of R band magnitude. The diamonds represent the extended 
objects, the crosses the stellar objects, and the squares the faint 
rriR > 21.5 identifications. The dotted lines show the upper and 
lower limits on colour for the expected range of spectral energy 
distributions for quasars. The solid lines indicate the expected 
colours for passively evolving ellipticals which are on the R band 
Hubble diagram found for radio bright GPS galaxies (Snellen et 
al. 1996a). 

than expected on the basis of the R band Hubble relation. 
Indeed the Hubble diagram is poorly determined for galax- 
ies with R > 23.0. The discrepancy of 1.5 magnitudes in 
R — K colour for a tur = 24.0 galaxy could be explained if 
it is not at a redshift of 1.3, but at a redshift of 1.0. This is 
well within the uncertainties of the Hubble diagram at such 
a high redshift. 

4.2 The Optical Morphologies and Environments 

It has been claimed by Stanghellini et al. (1993) that 75% 
(8 out of 11) of the GPS galaxies in the radio bright sam- 
ple with rriR < 20.5 have distorted isophotes, double nuclei 
or close neighbours. There are four galaxies in our sample 
which have tur < 20.5 and can be used to check this claim. 
The angular morphologies of B0830+5813, B1622+6630, 



B1819+6707 and B1946+7048 have ellipticities of 0.173, 
0.026, 0.070, and 0.031, with position angles of 29°, 61°, 
53°, and 51° respectively. Contour plots of the two bright- 
est galaxies, B0830+5813 and B1946+7048, seem to indicate 
that these are "relaxed" systems (as indicated by the regular 
symmetric contours). 

To check this an r 1 ^ 4 law was fitted to the galaxy profiles 
(figure^). The profile of B1946+7048 is well fitted by an r 1/4 
law, with an effective radius of R e ff = 13", but the profile 
of B0830+5813 clearly flattens in the inner 10 pixels. In this 
case an r 1 ^ law was fitted to the outer part of the profile 
{R e ff = 14"). This indicates that about half of the expected 
light is missing in the centre of the galaxy, possibly due to 
the presence of dust. Unfortunately the K band data on 
this galaxy is of too poor a quality to measure the colour 
gradient along the profile. Assuming that the light profile 
without dust followed the J? 1//4 law, then an estimate of the 
dust mass Mdust can be made. Following van Dokkum and 
Franx (1995), the dust mass corresponds to 

M dust = EAyTy 1 (6) 

where E is the area covered by the dust, Av is the average 
extinction in this area, and Vv ~ 6 x 10 -6 mag kpc 2 Mq 
visual mass absorption coefficient. The average extinction 
is estimated to be 0.8 mag within the radius Rdust = 1.4". 
The redshift of the galaxy is z=0.093 (Snellen 1997, Snellen 
et al. in prep) , hence this corresponds to an area of 30 kpc 2 
(Ho = 50 km/sec/Mpc, qo = 0.5). This gives a dust mass 
estimate of Mdust ~ 10 6 ' 5 Mq. This is comparable to dust- 
masses estimated in 3C galaxies (de Koff et al., 1998). 

The exposure times on B1622+6630 and B1819+6707 
were not long enough to construct reliable profiles such as 
for B0830+5813 and B1946+7048, but there are indications 
from the noisy profiles that they are also relaxed, although a 
nuclear point source might well be present in B1622+6630. 

Of the four brightest galaxies, two seem to have close 
neighbours but none of them show evidence for an ongoing 
galaxy interaction. There are two objects located at a dis- 
tance of ~ 30" from the core of B0830+5813, clearly within 
the envelope of the galaxy. Although these objects may well 
not be associated physically with the main galaxy and the 
result of a projection effect. B0830+5813 seems to be located 
in a cluster environment (fig. [|). There are 13 galaxies on 
the CCD image with magnitudes between the magnitude of 
the third brightest galaxy and 2 magnitudes fainter. This 
corresponds to an Abell richness class (Abell, 1958) of 2 — 3, 
assuming a modified Hubble profile (Sarazin, 1986) for the 
projected surface density of galaxies with a core radius of 
500 Kpc. The galaxies B1622+6630 and B1819+6707 seem 
to have close neighbours at a distance of respectively 2.0" 
and 1.5", but B1946+7048 has no neighbouring galaxies. 

4.3 The Optical Identifications and their Radio 
Spectra 

If the orientation unification scheme for extragalactic radio 
sources is correct, then it is expected that the redshift dis- 
tributions of radio galaxies and quasars are more or less the 
same. If differences occur then these must be due to lumi- 
nosity effects, for example, a change of quasar opening angle 
with radio power (eg. Urry and Padovani, 1995). The ob- 
served redshift distributions in radio-bright samples of GPS 



© 0000 RAS, MNRAS 000, 000-000 



Optical and Near-infrared Imaging of faint GPS Sources 7 




Figure 4. The angular profiles of B0830+5813 and B1946+7048, fitted with an i? 1 / 4 law. Note that the profile of B0830+5813 clearly 
flattens in the inner 10 pixels indicating that light is missing in the centre of the galaxy, which may be due to the presence of dust. 
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Figure 7. The observed radio peak frequency as function of 
R magnitude for the faint and/or extended objects, which are 
assumed to be the GPS galaxies. There is no correlation found. 



Figure 6. The radio peak frequency distribution for the faint or 
extended objects, assumed to be galaxies, and the bright stellar 
objects, assumed to be quasars. The two distributions are differ- 
ent at a 99% confidence level (KS-test). 



sources identified with galaxies and quasars are found to be 
very different, with the GPS quasars located at high red- 
shifts (2 < z < 4, O'Dea et al, 1991). However, two se- 
lection effects influence this result. Firstly, the GPS galax- 
ies are biased towards low redshift because spectroscopy is 
easier on low redshift galaxies. Secondly, the GPS quasars 
may be biased towards high redshift because their intrin- 
sic peak frequencies seem to be higher (eg. de Vries et al 
1997, Stanghellini et al 1996). If that is true, low redshift 
GPS quasars may fall outside the bright samples because 



their observed peak frequencies are too high. We have found 
evidence that the redshift distribution of GPS galaxies in 
the faint sample is indeed different from that of the GPS 
quasars (Snellen 1997). The faint sample is more suitable 
for this, because sources are selected on the optically thick 
part of their radio spectrum, and moving sources within the 
sample to lower redshift, thereby shifting the observed peak 
frequency to lower frequencies, will not make them disappear 
from the sample. The difference in redshift distributions is 
unlikely to be caused by a radio power effect, because it 
occurs both in the faint and bright samples. This seems to 
indicate that GPS galaxies and quasars are not unified by 
orientation and form two different classes of objects which 
happen to have the same radio spectral properties. 

As mentioned above, it has been found (Stanghellini et 
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Figure 5. The GPS galaxy B0830+5813, which is clearly located in a cluster environment 



al., 1996, de Vries etal. 1997) that the intrinsic peak fre- 
quency distributions of GPS quasars and galaxies are dif- 
ferent. However, the selection bias discussed above may be 
responsible for this. We conclude that it would be better 
to compare the observed peak frequency distributions of the 
galaxies and quasars. Although radio selection effects are se- 
vere in the observed frame, they are independent of optical 
host. If all optical counterparts of the faint GPS sources that 
are classified as "extended" or "faint" can be assumed to be 
galaxies, and all unresolved star-like identifications quasars, 
figure [] shows that there is a clear difference in the peak 
frequency distribution for these two groups. The hypothesis 
that the distributions of the quasars and galaxies are drawn 
from the same parent distribution can be rejected at > 99% 
confidence level using the Kolmogorov-Smirnov test. 

De Vries et al. (1997) have found an apparent absence 
of sources at high redshift with low peak frequencies, which 
they interpreted as a correlation between redshift or radio 



luminosity with peak frequency. We believe this is not a 
dependence of redshift, but a dependence on optical host 
because the quasars tend to be at higher redshift than the 
galaxies. To check this we investigated whether the observed 
peak frequency is a function of redshift for the GPS galaxies 
only. The redshifts are not known for the majority of the 
GPS galaxies, but the R band magnitude of GPS galaxies 
is a good indicator of the redshift, at least for the radio 
bright samples (Snellen et al. 1996a). Plotting the R band 
magnitude of the 'candidate' galaxies versus their radio peak 
frequency however, shows no correlation between these two 
quantities (fig. [?]). 



5 CONCLUSIONS 

A sample of 47 faint GPS sources selected from the West- 
erbork Northern Sky Survey has been imaged in the opti- 
cal and near- infrared, resulting in an identification fraction 
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of 87% with a reliability of 97.6% and a completeness of 
> 99.9%. The optical to near infrared colours of the faint 
and the extended optical identifications are found to be in 
agreement with passively evolving elliptical models, assum- 
ing that their redshifts can be determined from the R band 
Hubble diagram as defined for radio bright GPS galaxies 
(Snellen et al. 1996a), and taking into account the large un- 
certainty in the Hubble diagram at high redshift. 

The GPS sources optically identified with bright stellar 
objects, presumably quasars, have in the observed frame, 
on average, higher radio peak frequencies than the GPS 
sources identified with faint and extended objects, presum- 
ably galaxies. This indicates that there is a difference be- 
tween the true peak frequency distributions of GPS galaxies 
and quasars, because selection effects are the same for both 
type of optical counterparts. Such a difference in intrinsic 
peak frequency distribution may cause a selection bias, re- 
sulting in a different redshift distribution for GPS galaxies 
and quasars. However, it has been shown by this faint sam- 
ple (Snellen 1997) that this is not the case, and that both 
the redshift distribution and intrinsic peak frequency distri- 
bution for GPS galaxies and quasars are different. 

This seems to indicate that GPS galaxies and quasars 
are not unified by orientation and form two different classes 
of objects which happen to have the same radio spectral 
properties. 
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APPENDIX A 

Comments on individual sources 

B0400+6042 No optical identification has been found to a 
limiting magnitude of R — 22.1 and K — 18.5. 
B0436+6152 No optical identification has been found to a 
limiting magnitude of R — 22.1 and K — 18.5. 
B0441+5757 A very bright star is located only ~ 15" away 
to the north-west. There may be a second object to the south 
east very close ~ 2" to the identification. 
B0539+6200 An optical identification is found in R band, 
just at the detection limit. However, the object is clearly 
detected in K band. 

B0552+6017 No optical identification has been found to a 
limiting magnitude of R = 23.8. 

B0759+6557 This object has faint detections in both R 
and K band. 

B0826+7045 Two possible identifications are found at 0.8 
and 1.2 arseconds distance from the radio position. 
B1600+7131 No optical identification has been found to a 
limiting magnitude of R = 24.1 and I = 21.5. 
B1620+6406 A very bright star is located only 2 arcmin- 
utes to the south west. 

B1622+6630 Two objects are found within an angular dis- 
tance of 2". The may be neighbouring galaxies. 
B1639+6711 No optical identification has been found to a 
limiting magnitude of R = 24.1 and I = 21.5. 
B1647+6225 Two objects are found within 3.5 arseconds 
from the radio position. The fainest one is probably the real 
identification. 

B1655+6446 No optical identification has been found to a 
limiting magnitude of R = 24.1 and I = 21.5. 
B1807+5959 A possible identification is found at a 3cr 
magnitude level. 

B1807+6742 No optical identification has been found to a 
limiting magnitude of R = 22.7 and I = 21.5. 
B1819+6707 Two objects are found within an angular dis- 
tance of 1.5". The may be neighbouring galaxies. 
B1843+6305 This source has no identification to R > 22.6 
but is clearly detected in I band. 

B1942+7214 The optical identification is very faint in R 
but is clearly visible in I and J band. 

B1954+6146 The optical identification is a Red point 
source, with very faint extended emission in R band. 
B1958+6158 Two objects are found within 2.0 arseconds 
from the radio position. 
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Table 2. The optical parameters of the candidate identifications. 
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Figure 8. THIS FIGURE IS AVAILABLE AT: http://www.ast. cam. ac.uk~snellen The optical and near-infrared images 



Note that the angular scale is different for the K band images. 
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